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General stability criterion of two-dimensional inviscid parallel flow 
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A more restrictively general stability criterion of two-dimensional inviscid parallel flow is obtained 



analytically. First, a sufficient criterion for stability is found as either — /ii < 



< or < 



u" 



u-u s ^ " "* " ^ u-u s 

in the flow, where U s is the velocity at inflection point, fii is the eigenvalue of Poincare's problem. 
Second, this criterion is generalized to barotropic geophysical flows in /3 plane. Based on the criteria, 
the flows are are divided into different categories of stable flows, which may simplify the further 
investigations. And the connections between present criteria and Arnol'd's nonlinear criteria are 
discussed. These results extend the former criteria obtained by Rayleigh, Tollmien and Fj0rtoft and 
would intrigue future research on the mechanism of hydrodynamic instability. 



PACS numbers: 47.20.-k, 47.20.Cq, 47.20.Ft, 47.15. Ki 

The stability due to shear in the flow is one of the 
fundamental and the most attracting problems in many 
fields, such as fluid dynamics, astrophysical fluid dynam- 
ics, oceanography, meteorology et al. The shear insta- 
bility has been intensively investigated, which is to the 
greatly helpful understanding of other instability mech- 
anisms in complex shear flows. For the inviscid parallel 
flow with horizontal velocity profile of U(y), the general 
way is to investigate the growth of linear disturbances 
by means of normal mode expansion, which leads to 
the famous Rayleigh's equation . Using this equation, 
Rayleigh first proved a necessary criterion for insta- 
bility, i.e., Inflection Point Theorem. Then, Fj0rtoft Q 
found a stronger necessary criterion for instability. These 
criteria are well known and have been applied to un- 
derstanding the mechanism of hydrodynamic instability 
H H @ . Unfortunately, both criteria are only necessary 
criteria for instability, except for some special cases of the 
symmetrical or monotone velocity profiles. Tollmien |(J 
gave a heuristic result that the criteria are also sufficiency 
for instability in these special cases. 

The stable criteria also provide a way to categorize the 
velocity profiles of the flows. According to Rayleigh's cri- 
terion, the flows are stable if U"(y) ^ 0, where U"(y) 
denotes d 2 U/dy 2 . And according to Fj0rtoft's criterion, 
there is another kind of stable flows if U"(U — U s ) > 0, 
where U s is the velocity at the inflection point U'J = 0. 
Then if U"(U - U s ) < 0, can the flow still be stable? 
Is there another kind of stable flows besides the above 
flows? To answer these questions, a more restrictive cri- 
terion is needed. And the criterion itself is important for 
both theoretic researches and real applications. The aim 
of this letter is to obtain such a stability criterion, and 
other instabilities may be understood via the investiga- 
tion here. 

For this purpose, Rayleigh's equation for an inviscid 
parallel flow is employed 0, IE S IE • For a parallel 



flow with mean velocity U(y), the streamfunction of the 
disturbance expands as a series of waves (normal modes) 
with real wavenumber k and complex frequency uj = uj r + 
iu>i, where u>i denotes the grow rate of the waves. The 
flow is unstable if and only if u>i > 0. We study the 
stability of the disturbances by investigating the growth 
rate of the waves, this method is known as normal mode 
method. The amplitude of waves, namely </>, satisfies 



U" 



(U-c) 



0, 



(1) 



where c = uj /k = c r +ici is the complex phase speed. The 
real part of complex phase speed c r = oj r /k is the wave 
phase speed. In fact, Rayleigh's equation is the vorticity 
equation of the disturbance jj, . This equation is to be 
solved subject to homogeneous boundary conditions 



= at y = a,b. 



(2) 



There are three main categories of boundaries: (i) en- 
closed channels with both a and b being finite, (ii) bound- 
ary layer with either a or 6 being infinite, and (iii) free 
shear flows with both a and b being infinite. 

It is obvious that the criterion for stability is u>i = 
(cj = 0), for that the complex conjugate quantities </>* 
and c* are also a physical solution of Eq.Q and Eq.©. 
Multiplying Eq.Q by the complex conjugate <fi* and in- 
tegrating over the domain a < y < b, we get the following 
equations 



I^H 2 + fc 2 ||0H 2 ) + 



U"(U-c r ) 
WU-cP 1 



<t>\\ 2 ]dy = (3) 



and 



U" 



U-c 



\4>\\ 2 dy = Q. 



(4) 
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Rayleigh used only Eq.(0J to prove his theorem. Fj0rtoft 
noted that Eq.© should also be satisfied, then he ob- 
tained his necessary criterion. To find a more sufficient 
criterion, we shall investigate the conditions for ci — 0. 
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Unlike the former investigations, we consider this prob- 
lem in a totally different way: if the velocity profile is 
stable (c, = 0), then the hypothesis Ci ^ should result 
in contradictions in some cases. Following this, some 
more restrictive criteria can be obtained. 

To find a stronger criterion, we need to estimate the 
ratio of ||0'|| 2 dy to J \\(f)\\ 2 dy. This is known as 
Poincare's problem: 



(5) 



where the eigenvalue \i is positive definition for any 
(j) =/= 0. The smallest eigenvalue value, namely can 
be estimated as /ii > (^z-) 2 , like Tollmien [|| did. 

Then using Poincare's relation Eq.©, a new stabil- 
ity criterion may be found: the flow is stable if —fii < 
ifzir < everywhere. 

To get this criterion, we introduce an auxiliary function 
fiv) = uLjj i where f(y) is finite at the inflection point. 
We will prove the criterion by two steps. At first, we 
prove proposition 1: if the velocity profile is subject to 
~Mi < f(y) < 0, then c r ± U s . 

Proof: Since —\x\ < f(y) < 0, then 



Mi < 



U" 



U"(U - U s ) < U"(U - U„) 



(0) 



u-u. {u-u 8 y - (u-u s ) 2 + c 2 

Substitution of c r = U s and Eq. © into Eq. © results in 



/ 1 1 2 



U"{U-U S ) 
WU-cV ' 



<t>f)dy>Q. (7) 



This contradicts Eq.©. So proposition 1 is proved. 

Then, we prove proposition 2: if —fii < f(y) < and 
c r 7^ U s , there must be cf = 0. 

Proof: If cf 7^ 0, then multiplying Eq.Q by (c r — 
Ut)/ci, where the arbitrary real constant Ut does not 
depend on y, and adding the result to Eq.©, it satisfies 



u"(u - Ut) 

WU-cP 



\dy = 0. (8) 



But the above Eq.© can not hold for some special Ut- 
For example, let U t = 2c r — U s , then there is (U — U S ){U — 
U t ) < \\U-cf, and 



U " {u - Ut) =m {u -Z s){u - Ut) >-^ 



\\u-c\\* 

This yields 

rb 



/ 1 1 2 



||17- C ||5 



2 , U"{U-U t ),... u3 



(9) 



\\U-cf 



]U\\ J }dy>0, (10) 



which also contradicts Eq.©. So proposition 2 is also 
proved. 

Using 'proposition 1: if — /i X < f(y) < then c r ^ U s ' 
and 'proposition 2: if —fii < f(y) < and c r =/= U s 
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FIG. 1: Vorticity profiles within the interval —1 < y < 1. 
Profile 2(^ = cos(7rj//2), dashed) is neutrally stable, while 
profile 1 (^ = cos(y), solid) and profile 3 = cos(2y), dash 
doted) are stable and unstable, respectively. 



then Ci = 0', we find a stability criterion. If the velocity 
profile satisfies —fii < jjz^j- < everywhere in the flow, 
it is stable. Moreover, the above proof is still valid for 
< f(y), which is equivalent to Fj0rtoft's criterion. Thus 
we have the following theorem. 

Theorem 1: If the velocity profile satisfies either — fi\ < 



u" 



< or < 



u" 



the flow is stable. 



u-u._ . u-u. 

This criterion is more restrictive than Fj0rtoft's cri- 
terion. As known from Fj0rtoft's criterion, the neces- 
sary condition for instability is that the base vorticity 
^ = —U' has a local maximal in the profile. Noting that 
U"/(U — U s ) ~ C"/Cs near the inflection point, where 
£ s is the vorticity at inflection point, it means that the 
base vorticity £ must be convex enough near the local 
maximum for instability, i.e., the vorticity should be con- 
centrated somewhere in the flow for instability. A sim- 
le example can be given by following Tollmicn's way 
As shown in Fig^ there are three vorticity profiles 
within the interval — 1 < y < 1, which have local max- 
imal at y = 0. Profile 2 (U = — 2sin(7ry/2)/7r) is neu- 
trally stable, while profile 1 (U = — sin(y)) and profile 3 
(U = — sin(2y)/2) are stable and unstable, respectively. 

Moreover, the stabile criterion for the parallel inviscid 
flows can be applied to the barotropic geophysical flows 
in (3 plane, like Kuo did This is a generalized stable 
criterion, we state it as a new theorem. 

Theorem 2: The flow is stable, if the velocity profile 



satisfies either 



-Mi < < or < u " 



-r in the flow, 



u-u, ^ " " x " ^ u-u e 
where U s is the velocity at the point U" = (3. 

The criteria proved above may shed light on the inves- 
tigation of vortex dynamics. Both Theorem 1 and Fig^ 
show that it is the vorticity profile rather than the veloc- 
ity profile that dominates the stability of the flow. This 
means that the distribution of vorticity dominates the 
shear instability in parallel inviscid flow, which is essen- 
tial to understanding the role of vorticity in fluid. So an 
unstable flow might be controlled just by adjusting the 
vorticity distribution according to present results. This 
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FIG. 2: Growth rate uji as an function of wavenumber k, (a) 
for Uj — tanh(l.ly), (b) for U2 = tanh(1.3i/), both within the 
interval — 1 < y < 1. 



is an very fascinating problem, but can not be discussed 
in detail here. 

To show the power of the criteria obtained above, we 
consider the stability of velocity profile U — tanh(ay) 
within the interval —1 < y < 1, where a is a con- 
stant. This velocity profile is an classical model of mix- 
ing layer, and has been investigated by many researchers 
(see 0,11,13 an d references therein). Since U"(U — U s ) = 
—2a 2 tanh 2 (ay)/ cosh 2 (ay) < for — 1 < y < 1, it might 
be unstable for any a according to both Rayleigh's and 
Fj0rtoft's criteria. But it can be derived from Theorem 
1 that the flow is stable for a 2 < tt 2 /8. For example, 
we choose a± = 1.1 and 0:2 = 1-3 for velocity profiles 
Ui(y) and U^iy)- The growth rate of the profiles can 
be obtained by Chebyshev spectral collocation method 
with 100 collocation points, as shown in Fig[21 It is 
obvious that c, = for U\ and Cj > for U2, which 
agrees well with the criteria obtained above. This is also 
a counterexample that Fj0rtoft's criterion is not sufficient 
for instability. So this new criterion for stability is more 
useful in real applications. 

The present stable criteria give a affirmative answer 
to the questions at the beginning, i.e., there are some 
stable flows if U"(U — U s ) < 0. Based on the former 
criteria, the velocity profiles can be categorized as fol- 
lows: (i) without inflection point (Reyleigh's criterion), 



(ii) f(y) > (Fj0rtoft's criterion), and (iii) fix < f(y) < 
(present criterion). Then the flow might be unstable only 
for f(y) < /ii and f(y) changing sign within the interval. 
However, if f(y) changes sign somewhere within the in- 
terval [a, b], then the flow is stable. For that f(y) chang- 
ing sign implies U' s " — but U'J" ^ 0, so U" does not 
change sign near the inflection point. Thus Ci must van- 
ish in Eq.Q, i.e., the flow is stable for f(y) changing sign 
within the interval. In this way, the flow might be un- 
stable only for f(y) < fii somewhere, which will intrigue 
further studies on this problem. In fact, there are still 
stable flows if fi\ < f(y) is violated. 

Recall the proof of theorem 1, it is found that the fol- 
lowing Rayleigh's quotient /(/) plays a key role in deter- 
mination the stability of the flows. 

rm . S'[H'\\ 2 + f(y)U\\ 2 ]dy nn 
lit )=mm — ; (11) 

* />ll 2 
Noting that the proof of theorem 1 is still valid in the case 
of /(/) > 0. We have such result: the flows are stable if 
1(f) > 0. Though this criterion is more restrictive than 
that in theorem 1, it is inconvenient for the real applica- 
tions due to unknown value of Rayleigh's quotient 1(f). 
Theorem 1 is more convenient for the real applications 
in different research fields. 

The idea of categorization the velocity profiles of the 
flows may simplify the investigation of stability problem. 
It can be seen from Rayleigh's equation Eq.QJ that the 
stability of profile U(y) is not only Galilean invariant, but 
also independent from the the magnitude of U (y) due to 
linearity. So the stability of U(y) is the same as that of 
AU(y) + B, where A and B are any arbitrary nonzero 
real numbers. As the value of U"(U — U s ) in Fj0rtoft's 
criterion is only Galilean invariant but not magnitude 
free, it satisfies only part of the Rayleigh's equation's 
properties. On the other hand the value of U" /(U — U s ) 
satisfies both conditions, this is the reason why the crite- 
ria in both Arnol'd's theorems and present theorems are 
the functions of U"/(U — U s ). Since the stability of in- 
viscid parallel flow depends only on the velocity profile's 
geometry shape, namely f(y), and the magnitude of the 
velocity profile can be free, then the instability of invis- 
cid parallel flow could be called "geometry shape insta- 
bility" of the velocity profile. As the above investigation 
shows that the inviscid shear instability is only associ- 
ated with the geometry of velocity profile. The concept 
of " geometry shape instability" would be help in further 
investigations. This distinguishes from the viscous insta- 
bility, which is also associated with the magnitude of the 
velocity profile. 

As mentioned above, we have investigated the stabil- 
ity of the flows via Rayleigh's equation, while Arnol'd 
considered the hydrodynamic stability in a totally differ- 
ent way. He studied the conservation law of the inviscid 
flow via Euler's equations and found two nonlinear sta- 
bility conditions by means of variational principle. So 
what is the relationship between the linear criteria and 
the nonlinear ones? 
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It is very interesting that the linear stability criteria 
match Arnol'd's nonlinear stability theorems very well. 
Applying Arnol'd's First Stability Theorem to parallel 
flow, the stable criterion is < C\ < (U - U s )/U" < 
C2 < 00 everywhere in the flow, where C\ and C2 are 
constants. This corresponds to Ffortoft's criterion for 
linear stability, and is well known [3, [L0| . Here we find 
that Theorem 1 proved above corresponds to Arnol'd's 
Second Stability Theorem, i.e., the stable criterion is < 
C\ < —{U- U s )/U" < C 2 < 00 everywhere in the flow. 
Given C± — Arnol'd's Second Stability Theorem is 

equivalent to Theorem 1. Moreover, the proofs here are 
similar to Arnol'd's variational principle method. For 
the arbitrary real number U u which is like a Lagrange 
multiplier in variational principle method, plays a key 
role in the proofs. So that the above Theorem 1 is similar 
to Arnol'd's theorems. 

Unfortunately, Arnol'd's nonlinear stability theorems, 
though quite useful in the geophysical flows [T(|, are sel- 
dom known by the scientists in other fields. The main 
reason is that the proofs of Arnol'd's theorems are very 
advanced and complex in mathematics for most general 
scientists in different fields to understand. Although 
Dowling [ljj suggested that Arnol'd's idea need to be 



added to the general fluid-dynamics curriculum, his sug- 
gestion has not been followed even 10 years later. Com- 
pare with Arnol'd's theorems, the theorems proved here 
are equivalent in some sense but much simpler and easier 
to understand, therefore it is more convenient to use our 
new results in applications. 

In summary, the general stability criteria are obtained 
for inviscid parallel flow. These results, which are equiv- 
alent to Arnol'd's nonlinear theorems, extend the for- 
mer theorems proved by Rayleigh, Tollmien and Fj0rtoft. 
Based on the criteria, the velocity profiles are divided into 
different categories, which may simplify the further inves- 
tigations. In general, these criteria would intrigue future 
research on the mechanism of hydrodynamic instability 
and to understand the mechanism of turbulence. And it 
also sheds light on the flow control and investigation of 
the vortex dynamics. 

The author thanks Prof. Sun D-J at USTC, Dr. Yue P- 
T at UBC (Canada) and two anonymous referees for their 
useful comments. This work was original from author's 
dream of understanding the mechanism of instability in 
the year 2000, when the author was a graduated student 
and learned the course of hydrodynamic stability by Prof. 
Yin X-Y at USTC. 
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